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Abstract 

Background: Human epidermal growth factor receptor 2 (HER2) has an important role in cancer aggressiveness 
and poor prognosis. HER2 has been used as a drug target for cancers. In particular, to effectively treat HER2-positive 
cancer, small molecule inhibitors were developed to target HER2 kinase. Knowing that curcumin has been used as 
food to inhibit cancer activity, this study evaluated the efficacy of natural curcumins and curcumin analogs as HER2 
inhibitors using in vitro and in silico studies. The curcumin analogs considered in this study composed of 4 groups 
classified by their core structure, (3-diketone, monoketone, pyrazole, and isoxazole. 

Results: In the present study, both computational and experimental studies were performed. The specificity of 
curcumin analogs selected from the docked results was examined against human breast cancer cell lines. The 
screened curcumin compounds were then subjected to molecular dynamics simulation study. By modifying 
curcumin analogs, we found that protein-ligand affinity increases. The benzene ring with a hydroxyl group could 
enhance affinity by forming hydrophobic interactions and the hydrogen bond with the hydrophobic pocket. Hydroxyl, 
carbonyl or methoxy group also formed hydrogen bonds with residues in the adenine pocket and sugar pocket of 
HER2-TK. These modifications could suggest the new drug design for potentially effective HER2-TK inhibitors. Two 
outstanding compounds, bisdemethylcurcumin (AS-KTC006) and 3,5-bis((E)-3,4-dimethoxystyryl)isoxazole (AS-KTC021 ), 
were well oriented in the binding pocket almost in the simulation time, 30 ns. This evidence confirmed the results of 
cell-based assays and the docking studies. They possessed more distinguished interactions than known HER2-TK inhibitors, 
considering them as a promising drug in the near future. 

Conclusions: The series of curcumin compounds were screened using a computational molecular docking and followed 
by human breast cancer cell lines assay. Both AS-KTC006 and AS-KTC021 could inhibit breast cancer cell lines though 
inhibiting of HER2-TK. The intermolecular interactions were confirmed by molecular dynamics simulation studies. This 
information would explore more understanding of curcuminoid structures and HER2-TK. 
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Background 

Human Epidermal Growth Factor Receptor 2 (HER2) is 
one of the tyrosine l<inase receptors in EGFR family, 
which includes EGFR/ErbBl, HER2/ErbB2, HER3/ErbB3 
and HER4/ErbB4 [1]. Since there is no natural ligand 
specific to HER2, HER2 tends to form heterodimer with 
other ligand-induced members [2]. After dimerization, 
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the complex can trigger downstream pathways such as 
Ras/Raf/MAPK and PI3K/AKT pathways to increase 
cell growth, cell survival and cell differentiation [3,4]. 
Considerable evidencesshowed that HER2 over expression 
was involved in many types of cancer such as breast, 
ovarian, gastric and prostate cancers [5]. Therefore, HER2 
is considered as a drug target for cancer therapy focusing 
on inhibiting HER2 to reduce tumor growth. 

At present, there are two main approaches used to 
inhibit HER2, namely; monoclonal antibodies such as 
Trastuzumab, and small molecule inhibitors such as 
Lapatinib [6] and SYR127063 (called SYR for short) [7] 
targetingon tyrosine kinase domain (HER2-TK). Although 
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Trastuzumab can downregulate HER expression, 
cardiotoxicity and drug resistance can be found in 
Trastuzumab- treated patients [8,9]. Moreover, side effects 
such as diarrhea, rash or nausea can be observed in 
Lapatinib treatment [10]. Hence, new inhibitors are urgentiy 
required for HER2-overexpressed cancer treatment. 

Recently, in 2011, the first HER2-TK structure 
complex with pyrimidine compound SYR was released 
(PDB access 3PP0), providing the new understanding 
of the kinase structure [7]. Unlike the active- or 
inactive-conformations of EGFR-TK, HER2-TK configur- 
ation was somewhat in the middle of these typical confor- 
mations. It was named "the active-like conformation", due 
to, the orientation of the helix-aC-out, the DFG-in 
and unformed secondary structure of the activation 
loop. The second crystal of HER2-TK complex with 
TAK-285 (PDB access 3RCD) adopted the similar 
conformation as mentioned above [11]. 

Curcumin (also known as diferuloylmethane) is generally 
found as the major compound in rhizomes of turmeric 
plants; Curcuma longa Linaeusas yellow residue. It has 
been used as spice and ingredients in folk medicinal 
remedies in many Asian countries. The curcumin and its 
three natural analogs, curcumin II (demethoxycurcumin), 
curcumin III (bisdemethoxycurcumin) and cyclocurcumin- 
possess the remarkable pharmacological effects for centur- 
ies, such as anti-inflammatory [12,13], antioxidant [14], 
anti-carcinogenesis [15-18]. Moreover, curcumins is safe 
to use in high dose with non-toxic report [19,20]. Despite 
many advantages of curcumins, the poor stability and 
bioavailability profiles of curcumins are questionable 
when it comes to directly using crude curcumin as the 
potent and selective cancer drug. Many researchers have 
been focusing on the developing the curcumin ana- 
logs to enhance their stability and bioavailability. In 
particular, the novel series of curcumin-analog com- 
pounds have been synthesized and studied their effect 
in various cell targets [21-26]. They possess several prop- 
erties, potent activity against parasite in Trypanosoma 
and Leishmania species [21], antimycobacterial activity 
[22], inhibiting nitric oxide production from Lps-activated 
microglial cells [25] and estrogenic properties [23,24,26]. 
Thus, in this paper, we aimed to investigate the effect 
of this set of curcumin analogs on the HER2-TK 
activity using both experimental and computational 
methods. 

Curcumin has been shown to inhibit cancer growth by 
means of inhibiting several tyrosine kinases including 
EGFR, HER2, MAPK, phosphorylase kinase, pp60c-src 
tyrosine kinase, protein kinase C, and protein kinase A 
[18,27-34]. Furthermore, the curcumins can inhibit various 
types of cancer including breast cancer cells [15,28] and 
also induce the internalization of HER2 from cell surface 
[35]. Recendy, curcumin analog cyclohexanone has shown 



to selectively inhibit tyrosine kinase domain of EGFR, 
in vitro, in vivo and in silico studies [36] which reveals an 
opportunity for direct binding between curcumins and 
tyrosine kinase domains of other EGFR family members. 
Furthermore, the in silico screening of the natural database 
against HER2 kinase showed that such curcumins 
could interact with kinase through benzene rings for 
hydrophobic interactions and carboxyl groups for 
hydrogen bond formation [37]. 

In this study, we investigated interactions between 
curcumin analogs and HER2-TK by using virtual 
screening based on molecular docking in order to 
find potential compounds against HER2-TK. The hit 
compounds have been validated by different inhibitions 
between two types of breast cancer cell-lines with both 
HER2-overexpression and HER2-non-overexpression. 
Such findings might be useful for further development 
of curcumins as a new HER2 inhibitor in the future. 

Methods 

Computational procedures 
The preparation of ligand 

The two dimensional (2D) structure of 143 curcuminoid 
analogs were collected from the previous studies [21-26] 
(Additional file 1: Table SI). The ionization states, 
tautomers, stereochemistries, and ring conformations of all 
curcuminoid structures were calculated and OPLS-2005 
force field was applied using LigPrep module in 
Schrodinger package. These structures were used as an 
initial material during computational docking procedure 
to study interactions with the binding site of the HER2 
tyrosine kinase domain. 

The preparation of protein 

The atomic coordinate of HER2 tyrosine kinase domain 
(HER2-TK) was obtained from the crystallographic 
structure, accession no. 3PP0 in Protein Data Bank 
(PDB) [7]. This structure contains asymmetric dimer of 
HER2-TK complex with selective inhibitor HER2-TK, 
pyrrolo[3,2-d]pyrimidine-based potent, SYR. In order to 
perform the docking calculations, only chain A was 
selected as the target template. Another chain of 
HER2-TK as well as the co-crystalized ligand(s) and 
crystal water molecules were removed. Hydrogen 
atoms were assigned and parameterized with Optimized 
Potential for Liquid Simulation version 2005 (OPLS-2005 
force field) using the protein preparation wizard, which 
continuously minimized the whole structure by the Impref 
module in the Schrodinger package. 

Docking procedure using Glide standard precision mode 
(SP mode) 

The structures of protein and ligands were prepared as 
previously described. The OPLS-2005 force field was 
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applied to both protein and ligands. The complexes of 
HER2-TK and each curcuminoid, including the co- 
crystal ligand were generated with molecular docking 
approach using Grid-based Ligand Docking with 
Energetics (Glide)with standard precision mode (SP mode) 
[38,39]. The grid map was generated in Receptor Grid 
Generation by setting the center of the grid map around 
the catalytic site. Self-docking between HER2-TK and SYR 
was performed to validate all parameters before being 
applied to the study of interactions between HER2 
and curcumins. 

Post-docking analysis 

In order to handle considerable number of docking 
results, the sub-groups of modified core structure of 
curcuminoids were classified. Top ranks docking score of 
each sub-group were selected to further test in cell-based 
assay. In addition, the poor scores of each curcumin 
sub-groups were also chosen to be the control set in 
breast cancer cell-line assay. 

Molecular Interaction and stability in binding pocket 

All simulation steps were performed using the SANDER 
module of the AMBER 12 package and AMBER FF03 
force-field parameters [40]. The partial atomic charges 
of ligand were computed by using AMl-BCC method 
as implemented in the Antechamber module of the 
AMBER package. Their atom types and missing force field 
parameters were assigned based on the general AMBER 
force field (GAFF). Each complex was immersed in an 
isomeric truncated-octahedron box of TIP3P water 
molecules (10 A from the solute surface) and neutralized 
by additional Cr anions. The system was then minimized 
with the five-step procedure (Additional file 1: Table S2). 
All steps included 5,000 steepest-descent minimization 
cycles and 5,000 conjugate-gradient minimization cycles 
with different restraints on the protein structure. For 
the first step, harmonic restraints with a force con- 
stant of 5 kcal/(mol • A^) were used to immobilize the 
heavy atoms of protein coordinates, excluding hydro- 
gen atoms, at the starting positions, while solvent 
molecules were allowed to relax the unfavorable con- 
tacts with other solvent and solute molecules. For the 
second, third and fourth steps, harmonic restraints 
with force constants of 5, 1 and 0.5 kcal/(mol • A'^), 
respectively, were used to restrain the backbone of 
the protein. In the last step, the entire system was 
minimized with no positional restraints. 

With weak positional restraints on the protein 
(force constant of 5 kcal/(mol • A^)), all systems were 
heated from 0 to 300 K during a 200 ps MD simulations. 
After removing the restraints from the protein, we 
equilibrated the system with constant volume and set 
the constant temperature at 300 K for 500 ps. Note that 



we observed the equilibrium of energy (potential, kinetic 
and total energy), temperature, pressure, volume, density 
and RMSD before moving on to the production runs. The 
production MD simulations were performed from 30 ns 
while maintaining constant pressure and temperature. 
With a collision frequency of 1 ps"^, the temperature in 
all simulations was controlled by Langevin dynamics. 
Using an isotropic position scaling algorithm with a 
relaxation time of 2 ps, the pressure in NPT simulations 
was maintained at an average pressure of 1 atm. The 
random number generator was reseeded [41] for every 
simulation. A cut-off of 10 A and the particle mesh 
Ewald method were employed with the default parameters 
to calculate long-range non bonded interactions. With the 
tolerance parameter of 10"^ A, SHAKE constraints 
[42] were used to eliminate bond-stretching freedom 
for all bonds involving hydrogen, thereby allowing 
the use of a 2 fs time step. To monitor the stabilities 
of all systems, the Ca root-mean-square deviations 
(RMSD) were calculated. The RMSD of binding resi- 
dues within 5 A of the inhibitor were examined. The 
ptraj modules in the AMBER software were used to 
calculate the hydrogen bond occupancy and hydrogen 
bond distance between inhibitors and proteins [43,44]. 
All MD simulations were calculated on 22-node Linux 
High Performance Computer Cluster with 32 cores of 
AMD 2.2 GHz. 

The energy calculations were done as implemented in 
the MMPBSA.py script in AmberTools. The MM-PBSA 
approach is an acceptable method to compute the free 
energies of binding of ligands to proteins or to estimate 
the absolute free energies of molecules [43,44]. One 
hundred frames from the last 5 ns of each 30 ns MD 
studies were selected for the analysis of ligand binding 
energies, sampled at 50 ps intervals. Binding free energy 
was estimated from each energy terms as following 
equations (equation 1-4), 

AGmmPBSA = ^Gcomplex-(AGprotein + ^Giigand) (1) 

AGmmpbsa = <AEmm> + <AGsoiv>-<TAS> (2) 
<AGsoiv> = <AGpb> + <AGnp> (3) 

<AEmm> = <AEvdw> + <AEeIe> (4) 

<AGmmpbsa> is referred to final calculated MM-PBSA 
binding energy. It is described by the difference of 
AGcompiex by the summation of AGpiotein and AGngand 
(1). The free energy of each molecular system is given by 
the expression in the terms of equation (2). <AEmm> 
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is the total molecular mechanics energy in the gas 
phase, <AGsoiv> is a correction term (solvation free energy) 
of each system surrounded by solvent, and <TAS> is the 
entropy. <AEmm> includes electrostatic <AEeie>. and van der 
Waals <AEvdw> energies, while <AGsoiv> is the sum of 
electrostatic solvation energy <AGpb>, and the non- 
electrostatic solvation component <AGnp> (3-4). The 
polar contribution is calculated using PB model, while 
the non-polar energy is estimated by solvent access- 
ible surface area (SASA). In this study, <TAS> term 
was excluded. 

Experimental procedures 
Proliferation and viability assay 

Cell proliferation and viability were measured by tetrazo- 
lium 3-(4, 5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium 
bromide (MTT) assay. The reaction was catalyzed by mito- 
chondrial succinate dehydrogenease and requires NADH, 
which must be supplied by living cells. SKBR3 and MCF7 
cell line were seeded in flat-bottomed 96-well tissue culture 
plates as 1 X 104 cells/well/ 100 i^L and cultured overnight. 
Pure curcumin extracts and its analogs with different 
concentrations were dissolved in 100 |iL 10% FBS-RPMI 
1640 medium and added into the cells and incubated for 
48 h. Then 100 |iL of medium was removed and 10 |iL of 
MTT dye (Sigma-Aldrich; USA) was added, followed 
by 4hourof incubation. Subsequently, the supernatant 
from each well was aspirated off, leaving the purple 
form azan crystals. Optical density was measured by 
an ELISA micro plate reader at 540 nm with a refer- 
ence wavelength of 630 nm. Percentage of cell sur- 
vival was calculated by the formula below. Each assay 
was done in triplicate and the standard deviations were 
calculated. 
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ICio 



t 











MCF7 
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Molecular Dynamics 
simulations 

Figure 1 The rational workflow of this study, starting from in silico 
screening based on molecular docking and validated with cell line 
assays (SKBR3:a breast cancer cell line which over-expresses the 
HER2 gene product, MCF-7: a breast cancer cell line that absence 
of HER2 protein overexpression) and deep interaction study by 
molecular dynamic simulation. 



% Cell survival = Absorbance of treated well x 100 

Absorbance of vehicle control well. 

All curcumin analog compounds were synthesized and 
published [21-26] by the laboratory of Prof. Dr. Apichart 
Suksamran, Department of Chemistry and Center of 
Excellence for Innovation in Chemistry, Faculty of Science, 
Ramkhamhaeng University. 

Results and discussion 

Selection of the curcumin analogs and structure analyses 

Considering to the core structures of all 143 curcuminoid 
compounds, the middle-linear seven carbon linkage 
between two phenyl rings can be classified into four 
sub-groups, p-diketone, monoketone, pyrazole (N-N 
heterocyclic), and isoxazole (N-O hetorocyclic). Schematic 
diagram of the workflow is shown in Figure 1. All 143 
curcumin analogs were docked against tyrosine kinase of 
HER2 by Glide SP docking. A few compounds of each 



sub-group in the docking top ranks were selected. All 24 
chosen compounds were classified into four groups based 
on their core structures. From Table 1, compounds 
AS-KTCOOl to AS-KTCOll were classified into p-diketone 
group. Monoketone composed of two compounds, 
AS-KTC012 and AS-KTC013. For AS-KTC014 to 
AS-KTC017 were categorized as pyrazole curcumin ana- 
logs. The last group, isoxazole curcumin analogs consisted 
of seven molecules, AS-KTC018 to AS-KTC024. The 
docked conformations of all curcumin compounds were 
well oriented in the ATP-binding pocket of HER2-TK 
(Figure 2C and 2D). One phenyl-end of curcuminoids 
compounds oriented well in deep hydrophobic pocket, 
while another phenyl-end exposed to the open gate. They 
could form interactions either with MetSOl (adenine 
region) or with CysSOS (hydrophobic pocket II), depending 
on each configuration of analog (Figure 2D). According to 
the revealed three dimensional structure of HER2-TK 
adopted the active-like conformation, the binding cave 



Yim-im et al. BMC Bioinformatics 2014, 15:261 
http://www.biomedcentra l.com/1 47 1 -2 1 05/1 5/26 1 



Page 5 of 13 




AS-KTC009 



-8.66 



[21,22] 
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AS-iaC018 



-8.04 



[22] 
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Stayed in the tunnel shape rather than the opened cave 
as presented in EGFR-TK. This pocket possesses approxi- 
mately the volume of 475 A^. As the binding pocket 
adopted the tunnel-like shape, the curcuminoid structures 
are also in the linear cylinder shape that can fit into this 
pocket quite well. The hydrophobic-I site merged with the 
phosphate binding region (Ser783, Arg784, Leu785, 
Leu769, Gly770, Ala771, Met774 and Phe864) as the deep 
semi-closed site of ATP-binding pocket. This pocket is 
occupied by the trifluoromethyl-phenoxy fragment of 
SYR, the co-crystal ligand (Figure 2A and 2B). 
Furthermore, one of the phenyl rings and its hydrophobic 
substituent groups of curcuminoids compounds were 
perfectly fit to this deep hydrophobic pocket. The 
functional groups of modified middle-linear seven carbon 
linkage of curcurmin analogs pointed to the DFG motif, 
especially the adjacent residue Thr862. In addition, 
another phenyl-end and its hydrophobic substituent 



groups of curcumin analogs could span around adenine 
region to hydrophobic-II site (Figure 2C and 2D). 
This observation was found as the common interaction 
among curcuminoid compounds, conformed to the 
pyridine and amine fragments of SYR pointed to the 
same residues (Figure 2 A and 2B). These occupancy and 
interactions were considered important to increase 
selectivity and affinity for HER2 inhibitors [45-47]. 
However, the curcumin analogs could not interact 
with MetSOl as any other EGFR/HER2-TK inhibitors, 
which contain 4-anilinoquinazoline, pyrrolopyrimidine, 
pyrrolotriazine and cyanoquinoline cores. The carbonyl 
group of the linkage chain, from |3-diketone or monoke- 
toncurcumins, could form a hydrogen bond directly with 
one of the key residues Thr862 (Thr862 OH — OC curcu- 
mins). This interaction differed from that of SYR127063 
and other HER2 inhibitors that was a water-mediated 
hydrogen bond between N3 of quinazoline and Thr862 
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Figure 2 The illustration of A) two dimensional (2D) structure and B) three dimensional (3D) structure of SYR from the x-ray structure 
3PP0, while panel C) and D) present 2D and 3D structures of curcumin analogs in the binding pocket of HER2-TK from the 
docking results. 



[48]. As the direct hydrogen bond formation at this 
position was crucial to enhance binding affinity, N3 of 
quinazoline was modified to nitrile group to form the 
direct hydrogen bond [48]. These occupancy and 
interactions mimicked those of a heterocycUc core of 
HER2 inhibitor in an adenine ring of ATP [49]. All 
24 selected compounds have been processed in further 
experiments, to investigate the bioavailability profiles 
of curcumin analogs on two types of human breast 
cancer cell lines. 

The bioavailability profiles of curcumin analogs on two 
types of human breast cancer cell lines 

The HER2-TK inhibitory assay was performed using a 
commercialized HER2-TK, but unfortunately, activity of 
the HER2-TK could not be detected (data not shown). 
Moreover curcumins were stable and soluble in very low 
pH which is not suitable for HER2-TK to stay active. 
Therefore, the cell-based assay was chosen instead of the 
purified protein based assay to investigate the inhibitory 
effect of curcumin analogs to HER2-TK. In order to 
identify the curcumin analogs against HER2-TK, two dif- 
ferent types of human breast cancer cell lines, MCF-7 



and SKBR3 were performed. The SKBR3 is a breast can- 
cer cell line which over-expresses the HER2 gene prod- 
uct while the MCF-7 is in absence of HER2 protein 
overexpression [50-54]. Therefore, such compounds 
should be more effective against SKBR3 than MCF-7. 
The lapatinib which was tested on both breast cancer 
cell lines showed the IC50 on SKBR3 lower than MCF-7 
[55]. Twenty-four selected compounds were determined 
IC50 on both breast cancer cell lines by MTT assay as 
shown in Table 2. Although most of curcumin ana- 
logs have very similar structures, each analog showed 
different activities on both cells. In the ^-diketone 
group from Table 1, only AS-KTC006 showed the 
effective inhibiting to HER2-positive cancer cell-line 
(SKBR3), vice versa, lack of inhibiting to HER2- 
negative cancer cell-line (MCF-7). The core structure 
of this group composes of p-diketone which is modi- 
fied substituent from natural curcumins. On the other 
hand none of curcumin analogs in monoketone and 
pyrazole groups showed selectively inhibiting activity 
among both cancer cell-lines. In the last group, AS- 
KTC021 also presented the outstanding suppress on 
SKBR3 but not MCF7 among isoxazole analogs. The 
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Table 2 The inhibitory activity profiles of curcumin 
analogs on MCF7 and SKBR3 cells 





SKBR3 


MCF7 




SKBR3 


MCF7 


AS-iaC 


ICso 


ICso 


AS-KTC 




ICso 




(MM) 


(^M) 






(MM) 


001 


8.3 ± 0.6 


41 .9 ± 12.3 


013 


> 100 


> 100 


002 


1 3.0 ± 1.8 


444 ± 1 0.8 


014 


10.8 ±5.5 


9.9 ±3.5 


003 


24.9 ± 2.3 


8 1.6 ±26.0 


015 


30.9 ± 45 


154 ±3.8 


004 


> 100 


794 ± 9.8 


016 


42.6 ±5.5 


33.8 ±5.8 


005 


> 100 


> 100 


017 


21.3 ±3.8 


14.3 ± 1.9 


006 


154 ±3.9 


> 100 


018 


33.8 ±7.2 


224 ±7.7 


007 


7.9 ±2.5 


1 7.5 ± 45 


019 


> 100 


38.9 ±9.1 


008 


7.9 ±2.5 


36.7 ± 5.8 


020 


25 ± 5.5 


44.1 ±9.0 


009 


8.2 ± 04 


22.1 ±0.1 


021 


16.9 ±34 


> 100 


010 


> 100 


> 100 


022 


> 100 


> 100 


Oil 


> 100 


> 100 


023 


> 100 


24.3 ± 8.7 


012 


9.9 ± 1 .0 


14.3± 1.6 


024 


> 100 


> 100 



ICso for AS-KTC006 and AS-KTC021 in SKBR3 were 
15.4 and 16.9 \iM, respectively, and IC50 for both com- 
pounds in MCF7 were higher than 100 |iM. Since both 
AS-KTC006 and AS-KTC021 were selected from dock- 
ing results, they could inhibit the breast cancer cells 
through blocking of HER2-TK activities. Therefore, 
the AS-KTC006 and AS-KTC021 were chosen for fur- 
ther investigation the interaction mechanisms by molecular 
dynamics simulations. 

Molecular Interaction, stability binding free energy via 
MM-PBSA 

The molecular dynamics (MD) simulations were per- 
formed to examine the molecular interaction of 
the bothAS-KTC006 and AS-KTC021 curcuminoids in 
the ATP-binding pocket of HER2-TK. As mention in the 
previous section, the x-ray crystal structure of HER2- 
TK, 3PP0 has been used as the model reference template 
of this study. The complexes of HER2-TK-AS-KTC006 
and HER2-TK-AS-KTC021 which were constructed 
using molecular docking procedure have been used as 
starting coordinates for MD calculations. The root 
mean square deviations (RMSD) of all systems (back- 
bone atoms, ligand atoms and binding site atoms) 
seemed to fit nicely in the binding pocket of HER2-TK 
(Figure 3, Additional file 1: Table S3). The structure of 
curcuminoid distinguishes from other known tyrosine 
kinase inhibitors, which generally containing either 
quinazoline or pyrrolopyrimidine based structures 
[56]. The structures of curcumin analogs adopted the 
long thread with two knots at each end, resulting in 
freely flexible structure in the tunnel-like binding 



pocket. Interestingly, the MD results revealed that the 
hydrogen bonds between H-N atom of MetSOl 
(located on adenine region) and N-atom at position 
Nil of SYR existed about 94.73% along the entire 30 
ns of simulation time. On the other hand, the curcu- 
min analogs possess the interesting hydrogen bond 
pairs in both AS-KTC006 and AS-KTC021. Rather than 
forming hydrogen bond with MetSOl, the (3-diketone cur- 
cumin analog (AS-KTC006) formed hydrogen bonds to 
Thr862, Cys805 and Asp863, in the binding pocket of 
HER2-TK. The formation with these three amino 
acids existed approximately 34.67%, 10.17% and 
9.17%, respectively along the entire MD simulations 
(Table 3). Distances chromatograms of each pairwise 
atom from H-bond analyses were shown in Figure 4. 
One of these three residues is DFG-motif (Asp863), 
and another is adjacent residue of DFG-motif 
(Thr862). For AS-KTC021 (isoxazole curcumin ana- 
log), the hydrogen boding with Cys805 and Thr862 
existed about 35.97% and 5.37%, respectively, along 
the entire MD simulations. As illustrated in Figure 2C 
and 2D, one of the phenyl-end of curcumin com- 
pounds were well oriented in the ATP-binding pocket 
of HER2-TK, while the other-end exposed to solvent 
giving a chance to interact with Cys805, while the linker 
in the middle of cucurmin structure could interact 
with Thr862. The binding residues of SYR-HER2-TK 
systems appear to be stable along the MD simulations 
(Figure 3A-1 and 3A-2). The system started to converge 
since 15 ns of simulations time (Figure 3A). In addition, 
the simulation systems of AS-KTC006-HER2TK and 
AS-KTC021-HER2TK seem to converge after 25 ns of 
simulations time (Figure 3). From the molecular 
dynamic results suggested that AS-KTC006 had better 
binding affinity with HER2-TK than that of AS-KTC021. 
Considering at the curcuminoids structures, AS-KTC006 
is more flexible than AS-KTC021. In particular, AS- 
KTC0006 possesses P-diketonemoiety on both sides 
providing the flexibility of the molecule and allowing 
O-atom of ketone to interact with Thr862. The bind- 
ing energy calculations were performed to further 
understand the interactions of each system. It is to be 
noted that protein-ligand entropy contributions were 
not included in these free energy values since the 
present MM-PBSA are typically time consuming and 
unreliable. Considering the intermolecular interaction 
of the ligands with HER2-TK in contribution terms 
(Table 4), non-polar contributions, the summation of 
AEvDw arid AGnp, are significant with all ligand(s)- 
HER2TK systems. Interestingly, ligand AS-KTC006 
shows the most favorable electrostatic interaction with 
the binding pocket of HER2TK. There are the agree- 
ment between MM-PBSA binding energy calculations, 
docking and MD simulations. The consistent observations 
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C) AS-KTC021. 



were presented in all three different computational reported [57,58]. These worl<s showed that the vdW term 

approaches, namely, focusing on O-atom of l<etone to could be a major factor of the ligand-protein interactions; 

interact with Thr862 of the system AS-KTC006-HER2TK hence, the deep hydrophobic pocket would be the se- 

complex. Recently, the computational model of anti- lectivity pocket of HER2-TK [57,58]. Rather than 

HER2 ligands, the analogous of 4-anilinoquinazoline were focusing on the 4-anilinoquinazoline core structure 



Table 3 Conclusion of H-bonds between compounds and tyrosine kinase of HER2 



System 


Donor 






Acceptor 




% occupied 


distance (A) 


SYR-HER2TK 


MetSOl 


N 


H 


SYR 


Nil 


94.73 


3.1 28 ±0.1 5 




SYR 


01 


HI 


Asp863 


0D1 


10.67 


2.761 ±0.17 


AS-KTC006-HER2TK 


Thr862 


OGl 


HGl 


AS-KTC006 


024 


34.57 


2.851 ±0.18 




Cys805 


N 


H 


AS-KTC006 


025 


10.17 


3.191 ±0.18 




AS-KTC006 


021 


H39 


Asp863 


0 


9.17 


2.758 ±0.16 


AS-KTC021-HER2TK 


Cys805 


N 


H 


AS-KTC021 


022 


35.97 


3.142 ±0.1 8 




Thr862 


OGl 


HGl 


AS-KTC021 


021 


5.37 


3.100 ±0.20 
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(the known HER2-TK inhibitors), we focused on the 
curcuminoid core structure in this study. Both 
selected curcumin analogs form the distinguish inter- 
action moiety from the known inhibitors of HER2-TK. 
Furthermore, they occupied well in the deep hydrophobic 
pocket ofHER2-TK. 

Conclusion 

In the present study, we screened a series of curcumin 
compounds using a computational molecular docking. 
Then, the bioavailability assay of curcumin analogs, 
were conducted on two types of human breast cancer cell 



lines to select the specific active HER2 kinase 
inhibitors. The results suggested that bisdemethylcurcumin 
compound (AS-KTC006, CAS no. 60831-46-1) and 3,5-bis 
((E)-3,4-dimethoxystyryl)isoxazole (AS-KTC021, CAS no. 
1118765-46-0) could inhibit breast cancer cell lines 
though HER2-TK. In addition, the intermolecular studies 
from MD simulation suggested that both selected curcu- 
min analogs form the distinguish interaction moiety 
from the known inhibitors of HER2-TK. MM-PBSA 
binding calculation suggested that non-polar contributions 
are not only significant with all ligand(s)-HER2TK systems 
but also a major factor of the ligand-protein interactions. 



Table 4 Individual terms of MM-PBSA binding energy (kcal mol entropy term excluded 



System 


<AEvdw> 




<AEe|e> 




<AGpb> 




<AG„p> 




<AGsoi„> 


<AGjvimPBSA' 


Nonpolar/ 
hydrophobic 


Polar/electrostatic 


SYR-HER2TK 


-67.85 


(0.32) 


-16.38 


(0.53) 


45.48 


(0.44) 


-43.74 


(0.13) 


1 .73 (046) 


-82.50 (0.45) 


-111.59 (0.45) 


29.10 


(0.97) 


AS-ia"C006-HER2TK 


-42.81 


(0.28) 


-57.70 


(046) 


7627 


(0.37) 


-34.02 


(0.08) 


42.25 (0.37) 


-58.26 (0.39) 


-76.83 (0.36) 


1857 


(0.83) 


AS-KTC021-HER2TK 


-51.24 


(0.30) 


-18.15 


(0.63) 


37.15 


(0.63) 


-34.75 


(0.16) 


241 ( 0.60) 


-66.98 (0.38) 


-85.99 (0.46) 


19.00 


(1.26) 



<AEvdw^ 3nd (AEeie> - van der Waals and electrostatic contributions to binding energy. 
<AGpb> and <AGnp> - electrostatic and nonpolar contributions to the solvation free energy. 
<AGmivipbsa> " fin^l calculated MM-PBSA binding energy. 
Nonpolar contribution = AEydw + AGnp," polar contribution = AEgie + AGpb. 
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Additional file 



Additional file 1: Table S1-S3. Table SI. List of 143 curcuminoids 
compounds, which used in this present study. Table S2. Explanation of 
each simulation steps of minimization and molecular dynamics 
simulations. Table S3. List of binding residues of each system, which 
mentioned in Figure 3. 
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Growth Eactor Receptor; HER2-TK: Tyrosine kinase domain of HER2; RMSD: 
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Glide: Grid-based Ligand Docking with Energetics; SP mode: Standard precision 
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